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WHAT IS CLAIMED IS: 

1 . A method for determining vibration amplitude limits to detect faults in 
mechanical equipment, comprising: 

estimating a data probability distribution based on data for the mechanical 
equipment; and 

utilizing the data probability distribution to calculate the vibration 
amplitude limits. 

2. The method of Claim 1 further comprising removing outlier data. 

3. The method of Claim 2 further comprising calculating the vibration 
amplitude limits as a function of frequency for a substantial portion of the frequency 
spectrum. 

4. The method of Claim 1 wherein the data probability distribution is 
calculated using statistics and historical data of the mechanical equipment. 

5. The method of Claim 4 further comprising specifying importance levels 
for certain frequencies. 

6. The method of Claim 5 wherein the certain frequencies comprise 
frequencies for at least one of a motor, a compressor, or a gear. 

7. The method of Claim 6 further comprising obtaining vibration spectra 
comprising individual spectrum for the mechanical equipment from a database. 

8. The method of Claim 7 further comprising calculating a frequency for the 
individual spectrum and identifying the individual spectrum having a smallest number of 
frequency lines. 
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9. The method of Claim 8 further comprising calculating noise bandwidths 
and a largest noise bandwidth. 

10. The method of Claim 9 further comprising collecting vibration data from 
all spectra in a given frequency range. 

1 1 . The method of Claim 4 wherein the data probability distribution is 
calculated using a kernel density method. 

12. The method of Claim 1 1 wherein the kernel density method comprises 
calculating conditional kernel density. 

13. The method of Claim 12 wherein calculating conditional kernel density 
comprises estimating an unknown probability density for a given dataset. 

14. The method of Claim 13 wherein the probability density estimate at a 
point x for a one-dimensional dataset with n data points is given by: 



1 n 

P( x )=—r L K 



^ X — Xj ^ 



where, xj is the / h observation of dataset X, h is a bandwidth that characterizes a spread of 
the kernel, and k (•) is a kernel density function that is symmetric and satisfies the 

condition: £ K(u)du = 1 . 

15. The method of Claim 14 wherein the kernel density estimate is a two- 
dimensional kernel density estimate utilizing frequency and amplitude directions of the 
frequency spectrum. 
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16. The method of Claim 15 wherein a rf-dimensional kernel density estimate 
is generally written as: 



where K (w) is a d-dimensional kernel, H is a bandwidth matrix, and | ■ | denotes a matrix 
determinant. 

17. The method of Claim 4 further comprising detecting one or more faults in 
the mechanical equipment. 

18. The method of Claim 1 wherein the mechanical equipment comprises one 
or more HVAC chillers. 

19. A method for detecting faults in a chiller based on vibration amplitude 
limits, comprising: 

calculating vibration amplitude limits of the chiller using statistics and 
historical data for the chiller; 



generated; 

wherein the vibration amplitude limits are calculated as a function of 
frequency for an entire frequency spectrum. 

20. The method of Claim 19 further comprising removing outlier data. 

2 1 . The method of Claim 20 wherein the at least two-dimensional density 
estimate utilizes frequency and amplitude directions of the frequency spectrum. 

22. The method of Claim 21 wherein the at least two-dimensional density 
estimate is a ^-dimensional kernel density estimate. 




estimating an at least two-dimensional density estimate; and 
weighting the historical data based on when the historical data was 



-24- 



Atty. Dkt. No.: 081445-0359 



23. The method of Claim 22 wherein the ^/-dimensional kernel density 
estimate for point x of a dataset with n data points is given by: 



where, xj is the j observation of the dataset, K (u) is a ^-dimensional kernel, H is a 
bandwidth matrix, and | • | denotes a matrix determinant. 

24. The method of Claim 22 further including obtaining vibration spectra 
comprising individual spectrum for the chiller from a database. 

25. The method of Claim 24 further comprising calculating a frequency for 
the individual spectrum and identifying an individual spectrum having the smallest 
number of frequency lines. 

26. The method of Claim 25 further comprising calculating noise bandwidths 
and a largest noise bandwidth. 

27. The method of Claim 26 further comprising collecting vibration data from 
all spectra in a given frequency range. 

28. The method of Claim 19 further comprising calculating a conditional 
kernel density. 

29. The method of Claim 28 wherein calculating the conditional kernel 
density comprises estimating an unknown probability density for a given dataset. 
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30. A method for determining vibration amplitude limits of a mechanical 
device comprising: 

identifying a mechanical device and a frequency range for a spectrum to 

be analyzed; 

retrieving vibration spectra comprising individual spectrum for the 
mechanical device and the frequency range; 

calculating frequency for the individual spectrum; 

identifying the individual spectrum with a smallest number of frequency 

lines; 

calculating noise bandwidths and a largest noise bandwidth; 

removing outlier data; 

calculating conditional kernel density; and 

calculating vibration amplitude limits to detect faults in the mechanical 

device. 

3 1 . The method of Claim 30 wherein the mechanical device comprises a 
chiller for an HVAC system. 

32. The method of Claim 30 wherein the vibration spectra for the mechanical 
device and the frequency range is obtained from a database. 

33. The method of Claim 32 wherein calculating conditional kernel density 
comprises estimating an unknown probability density for a given dataset. 
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34. The method of Claim 33 wherein the probability density estimate at a 
point x for a one-dimensional dataset with n data points is given by: 



where, xj is the j observation of the dataset, h is a bandwidth that characterizes a spread 
of the kernel, and k (•) is a kernel density function that is symmetric and satisfies the 

condition: \ K(u)du = 1 . 



35. The method of Claim 33 wherein the kernel density estimate is at least a 
two-dimensional kernel density estimate utilizing frequency and amplitude directions of 
the frequency spectrum. 

36. The method of Claim 35 wherein a ^-dimensional kernel density estimate 
is given by: 



where K (w) is a ^-dimensional kernel, H is a bandwidth matrix, and | • | denotes a matrix 
determinant. 





-27- 



